Shot noise suppression at ID hopping 
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We have carried out a preliminary analysis of shot noise at 
hopping, focusing on uniform ID arrays of sites separated by 

tunnel barriers. The results show that at low temperatures 
the low- frequency density of the shot noise varies from to 
1 of the Schottky value, depending on the geometry, electron 
density, and Coulomb interaction strength. An interesting 
feature is dependence of the current spectral density at 

intermediate frequencies, which reflects self-similarity of the 
fluctuations at different size scales. 



I. INTRODUCTION 

Nonequilibrium fluctuations in mesoscopic systems can 
present additional information that is not reflected in 
their dc transport characteristics. This is one of the rea- 
sons why "shot noise" (i.e. nonequilibrium fluctuations 
of current with constant or nearly-constant spectral den- 
sity at low frequencies) has attracted so much attent 
in mesoscopics during the last decade - see, e.g., Refs.u 

An additional motivation Jor the present paper was 
provided by the observationa that the smallness of the 
shot noise is a necessary condition for quasi-continuous 
electron transfer. More exactly, for an external observer a 
conductor provides effectively Ohmic (quasi-continuous) 
conduction only if the so-called Fano factor 



F = SiiO)/2e{I) 



(1) 



(where Sj (0) is the low- frequency density of current fluc- 
tuations, and (/) is the average current) is much lower 
than 1. If simultaneously the resistance of such a sample 
is sufficiently high, and its stray capacitance is low. 



R > h/e^, C < eVr, 



(2) 



it may be used for resistive coupling in single-electron de- 
vices. Since using resistively-coupled devices is one of the 
very few options available to avoid the forbidding prob- 
lem of random |-hackground charge in single-electronics 
(see, e.g., RefsJj'tl), the search for systems with quasi- 
continuous conduction is important for possible future 
applications of single-electron devices in integrated cir- 
cuits. 

Shot noise has been extensiv|ely analyzed fcp-|rftetallic 
conduction (in both ballisticLTEl ancL.diffusivell3nlJ lim- 
its) and for single-electron tunncling.li3 Il3 Unfortunately, 
metallic conductors can satisfy the condition F <C 1 only 



if they are much_,kmger than the electron-phonon in- 
teraction length.EjO As a result Eqs. (||) can be prac- 
tically met|_Dnly at very low temperatures - see, e.g., 
experimentEj. The same is true for single-electron cir- 
cuits (UkjH! ID or 2D arrays) with their relatively large 
islandsl'Q. 

Much higher resistance R at small sample length (and 
hence small jC*) is typical for hopping conductors - see, 
e.g., Refs.E]^El. Naively, one might think that since the 
hopping transport is due to discrete single electron tun- 
neling events ("hops"), the shot noise should be close to 
the Schottky value (S'/(0) = 2e(/), i.e., F = 1). How- 
ever, this argument is obviously not true, since it could 
also be applied to a ID series of N tunnel junc- 

tions. A simple "circuit" theorytSo (see also Appendix) 
shows that for such an array the Fano factor can be very 
small: 



F - 1/iV < 1. 



(3) 



The physical reason for this fact is that the noise orig- 
inating from each junction is strongly shunted by the 
junction resistance, which is much smaller than the total 
resistance of other junctions. 

Thus, there is hope of having the shot noise at hop- 
ping suppressed well below the Schottky value as well. 
However, the real picture of hopping is complex, and the 
noise may be much higher than the simple estimate given 
above. For example, mutual correlation of the hopping 
events, exponentially broad distribution of their rates due 
to sample randomness, and the pesfjolative character of 
transport paths in 2D and 3D caseaH3^Eil may all be impor- 
tant factors. Until recently, the situation was virtually 
unexplored: the few publications on the theory of noise 
in hopping we are aware of (see, e.g., Ref.L3 and refer- 
ences therein) concentrate-pn 1// noise rather than on 
the broadband shot noise. E3 We are also unaware of any 
experimental studies of noise at hopping at frequencies 
high enough to avoid 1// noise dominance. 

The goal of this paper is to develop an initial picture of 
shot noise at hopping. We will focus on the ID case, and 
assume uniformity of hopping conditions between all the 
sites. (A brief analysis of nonuniform systems and higher 
dimensions is given in Discussion.) In principle, ID hop- 
ping may be implemented experimentally using a linear 
array of quantum dots between two external electrodes 
(Fig. I^a). Besides this geometry, we will also consider a 
somewhat artificial model of hopping on a ring (Fig. |l|b), 
at least because problems with periodic boundary condi- 
tions are traditional in theoretical studies of hopping. Be- 
sides that, since such models automatically conserve the 
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total electron number, they may crudely mimic "open" 
models (Fig. |^a) with considerable Coulomb interaction 
without its explicit account. 

Throughout our analysis we will assume that the elec- 
tron states localized at each site are non-degenerate, so 
that each site may be occupied with just one electron, or 
none. This model can be viewed as a special case, of the 
"orthodox" theory of single-electron tunnelingtfl when 
the background charge of each island is close to — e/2, so 
that energies of two charge states (n = and n = 1) are 
close to each other while other states are far beyond the 
available energy range. So, the well-ideyeloped theory p£ 
noise based either on Fokker-Planckli3~t3 or LangevinE2l 
approach can be directly applied to any hopping struc- 
ture with arbitrary electron-electron interaction. How- 
ever, these approaches involve taking into account an ex- 
ponentially large number of charge configurations, thus 
limiting practical calculations to relatively small struc- 
tures, N ^ 20. This is why for the numerical results 
we have used the Monte-Carlo approach, sipiilar to that 
used for simulations of transportEJ and noiseEju at single- 
electron tunneling, with the corresponding restriction of 
the site state number. 

It is instructive to compare the results for the shot 
noise at ID hopping and at tunneling in ID array of tun- 
nel junctions. Some formulas necessary for this compar- 
ison are derived in Appendix for the cases of negligible 
and small charge discreteness effects. 



where the factors \i depend on the structure geometry 
and can be expressed via its electrostatic matrix. (In 
general, these coefficients are different for the left and 
right electrodes.) In the simplest case of a ID array be- 
tween two infinite parallel metallic plates, — a.i/L, 
L = X^j^ij where a; is transport direction component 
of ith hop vector. In this work, we will use this formula, 
with Gi = L/N = const (i.e., Xi = 1/-/V), even for the ring 
geometry (Fig. ^), though this model does not have any 
electrodes. This assumption is not critical for the Fano 
factor which does not depend on Ai , since at low frequen- 
cies the spectral densities of all currents li and / coincide. 
(A simple proof of this statement may be obtained from 
the spectral density definition 



Si{lu) = lim -(I r I{t)e'^'dt\^) 

T^oo T Jo 



(6) 



in the limit ^ 0, using the condition that the charge 
cannot accumulate indefinitely inside the array). In the 
opposite limit of high frequencies (much higher than the 
average tunneling rate, though still much lower than the 
reciprocal "time of tunneling" , which is considered in- 
finitely short in our theory) , the spectral densities of cur- 
rents li and / are typically different, and obey a simple 
formula. In fact, in the high-frequency limit all tunnel 
events are effectively uncorrelated and the phases of fac- 
tors e'}cp{iLut^) in Eq. (^ are random. From this, we 
obtain 



II. SOME GENERAL RELATIONS 



5z,(^) = 2e((/+) + (/r)), Si{^) 



MSuioo). (7) 



In the hopping limit, where quantum interference be- 
tween states before and after each hopping event is 
neglected i beca use of the inelastic nature of electron 
transport ,E2lEil site occupation numbers may be consid- 
ered as random classical variables. If we are not inter- 
ested in extremely high frequencies (when the finite pho- 
ton energy becomes important), current Ii{t) flowing be- 
tween the {i — l)th and ith site may be considered as a 
sum of infinitely short pulses: 



m = I+{t)-lr{t), 



/±(t)=^e<5(t-i±), (4) 



where {t^ ) is the time of kth hop in the positive (neg- 
ative) direction between the sites. In the "open bound- 
ary" problem (Fig. ^), with a fixed voltage across the 
sample, we may alsp-jConsider currents I{t) flowing in 
external electrodesD'Eil These currents contain contribu- 
tions not only from the hops to and from the electrodes, 
but also the polarization charge changes (displacement 
current contribution) due to hops between internal sites: 



N 



(5) 



It is easy to see that for the current through one bar- 
rier 5/. (cx))/2e(/) > 1, while for the external current 
Si{oo)/2e{I) > 1/N. We will mostly be interested in the 
readily measurable quantity Si{uj) and its low frequency 
value Si{0). 

For the numerical (Monte-Carlo) calculations of the 
spectral density we have directly usedE£l Eq. (^). The 
time period t is chosen to be sufficiently long and the av- 
eraging is done over many such time periods. In practical 
calculations, it is important to keep the product u!t/2tt 
integer in order to avoid numerical inaccuracy at low fre- 
quencies, and it is convenient to calculate simultaneously 
the spectral density at several overtones of certain basic 
(low) frequency. For several figures we have also used 
the newly developed method for the calculation of spec- 
tral density, which gives much faster convergence; this 
method will be described elsewhere. 



III. CIRCULAR ARRAY 

A. The model 

We start with the auxiliary problem of hopping of 
a fixed number (M) of electrons on a uniform ring of 
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N > M sites. The electron may hop to either of the 
neighboring sites, i.e. either clockwise (with a probabil- 
ity rate of r+) or counterclockwise (with rate < r+), 
but only if the accepting site is empty. The rates 
should satisfy the Gibbs relation 



r-/r+ = exp(-iy/T), 



(8) 



where W is the energy difference between the neighbor- 
ing sites. (Due to the circular geometry, a conceptually 
sound, though impractical, way to create this difference 
is to increase the magnetic flux through the ring area at 
a constant rate. However, we consider the circular array 
mostly as a simplification of the realistic linear array.) 

For the comparison of the current noise with the 
Nyquist formula we will need the total resistance of the 
ring which is naturally defined as 



i?s = V/{I), 



(9) 



where the total "voltage" V is defined as NW/e (the 
dependence of the tunneling rate on W can be arbitrary). 

In the final part of our analysis we will include the par- 
ticle interaction following the unscreened Coulomb law, 
so that the potential energy of the system is 



[/{ri, r2, .., .rM}= ^ -p— ^ 



r,: - r,- 



(10) 



where rii is the site occupied by the ith electron. The 
interaction is included into the model by adding the cor- 
responding change of [/ at a hop to that (±VF) describing 
the external field. Since in this case the tunneling rates 
are no longer constant, we will need to specify an explicit 
relation (W) . In this case we will assume 



F = F+ - F- = W/eRo, 



(11) 



where Rq gives the scale of the effective resistance of a 
tunnel barrier between adjacent sites. 



B. Single particle limit 

Let us assume Then the current I{t) con- 

sists of uncorrelated pulses, each transferring the charge 
±e/-/V, with rates F+ and F~, respectively. This is equiv- 
alent to the conventional case of one tunnel junction with 
the electron charge substituted by e/N, hence 

(/)=e(F+-F-)/7V, (12) 

and the spectral density is frequency-independent ,0 
Si (to) = SiiO), with 

= ^ (r+ + r~) = 2e(/) ^ coth ^ . (13) 



Figure g shows the corresponding Fano factor F = 
coth(W^/2T), as a function of W. In thermodynamic 
equilibrium, W ~ 0, the noise satisfies the Nyquist for- 
mula. Si ~ 4T/Rj2, which remains valid while W ^T. 
At W ^ the Fano factor tends to infinity because 
the average current vanishes while the equilibrium ther- 
modynamic fluctuations still produce a finite current 
noise. For 1 the Fano factor crosses unity at 

~ VFc = 2r/A^ < T. Let us emphasize that since 
Wc ^ T, the noise at this crossover is still due to ther- 
modynamically equilibrium fluctuations. Finally, if the 
applied fleld is high {W ^ T), the Fano factor is low: 



F = 1/N. 



(14) 



Thus, as a matter of principle the shot noise suppression 
at hopping may be really very strong (proportional to 
the array length, just as in tunnel junction arrays). Now 
let us examine how this suppression is affected by various 
factors. 



C. Low temperature, no interaction 

At T < ly (i.e., F « F+ > F") and in the absence 
of Coulomb interactions (e^/a ^ W), but for arbitrary 
electron density p = M/N our model is reduced to the 
so-called Asymmetric Simple Exclusion Process (ASEP) 
model which has been extensively studied during the past 
few years - for a review see RefrZl. Within this model, all 
N\/Ml{N — M)l possible charge conflgurations of the sy^ 
tem have equal probability for the arbitrary N and M.EZl 
From this fact, the average current is readily calculated 
to equal 



(I) 



M N -M 
TV - 1 ' 



so that for a large system {N, M oo) 
(/) =eFp(l-p). 



(15) 



(16) 



Notice that these expressions (as well as those below) are 
obviously symmetric over the transformation p <-5- 1 — p 
which interchanges electrons and holes. From Eq. (p^), 
the maximum value {I) max = eF/4 of dc current is 
achieved at p = 1/2, which is a trade-off between in- 
creasing concentration p and decreasing average velocity 
F(l — p) of each electron (in hops per unit time) because 
of other electrons blocking its hops. 

Equation ( |l6|) is exactly the result which could be an- 
ticipated in the complete absence of correlation between 
the hops. However, in fact these correlations do exist, 
as revealed, for example, by the spectral density of the 
current. Figure |^ shows the result of numerical calcu- 
lation of Si{uj) using the Monte-Carlo approach for two 
concentrations, p — 0.3 and p = 0.5, and several val- 
ues of the array length N. The frequency dependence of 
the spectral density is obviously not flat as it would be 
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in an uncorrelated case. With increasing N the spectral 
density decreases and forms three distinct regions as a 
function of frequency: low and high frequency saturation 
regions and almost power-law decay in between. 
At high frequencies, in accordance with Eq. ^ 



S'/,(cx3) = 2e(/), Si (oo) ^2e{I)/N, 



(17) 



the suppression of the external current fluctuations is 
maximal. Notice that the frequency ujh of the crossover 
to this limit apparently does not depend on N , while the 
low-frequency crossover occurs at frequency uji which de- 
creases with N crudely as uji cx N^^^'^. (In tunnel junc- 
tion arrays, uji scales as N~'^ - see Appendix). rlJie zero 
frequency limit has been followed analyticallyEII giving 
the following Fano factor: 



F = 



rl/2 



N 



1/2 



(an analytical formula is also availablecZl for arbitrary 
N and M). Figure || shows that at large TV the fre- 
quency dependence of the current spectral density in 
the intermediate frequency range approaches power law: 

This dependence may be interpretedj-as a consequence 
of the self-similarity of the fluctuationa^J which occur at 
any site number scale L, within the interval 1 L <^ N . 
In order to explain the w^^/'^ scaling, let us assume that 
Eq. (^6|) is applicable to long-wave density perturbations 
in our system and introduce two velocities (measured in 
sites per second) of their propagation. _ 

The first of them, the sound ("shock"LZl) velocity 



To calculate fluctuations I{t) at a frequency w <C F, 
we can integrate Eq. (^6|) over the whole circle taking 
into account local density fluctuations Sp. Since we have 
assumed uniform Xi in Eq. (p|) and the total number of 
electrons does not fluctuate, jp{x)dx ~ AI, the contribu- 
tion from the linear term SI = eVSp vanishes. However, 
the current fluctuations do appear in the next, quadratic 
term of Eq. ([l^): 51 = — er((5p)^, which describes the 
"rectification" of density fluctuations. 

The density fluctuations at the size scale L {1 ^ 
L <C N) are described by the binomial distribution, 
giving the variance {{5p)'^) = p{l — p)/L. Hence, the 
typical relaxation bandwidth of these fluctuations (in 
the frame rotating with velocity Vs) is ujl — \^v\/L ~ 
F[p(l — p)]^/^L~'^/^, and the corresponding spectral den- 
sity is 5p(wi) ~ {5pf/LUL ^ [p(l -p)L]i/2/r. According 
to the standard theory of noise rectification (see, e.g.J^j), 
for7V,M^oo (18) Si{ujl) can be estimated as {N/L)e'^T^[Sp{ujL)?ujL, 

where the first factor accounts for N/L virtually inde- 
pendent fluctuating regions. Combining these estimates 
and eliminating L (as a function of ujl), we finally obtain 



SjiLu) ^ (a;/27rr)-V3 
MTi N [^(^ 



P)] 



2/3 



(21) 



V, = (1 - 2p)r, 



(19) 



(19) 



can be found from the obvious continuity equation 
d{ep)/dt — —dl/dx, where / and p are understood 
in the sense of "local" averages over 1 ^ SN <C N 
sites and x is the site number considered as a contin- 
uous coordinate. Since these averages are related by Eq. 
(p^), for small deviations from equilibrium we get dp/dt 
= —(1 — 2p)T dp/dx^ i.e. an equation describing linear 
waves propagating with the speed given by Eq. 
Notice that the sound velocity vanishes at half-fil 
p = 1/2, and is negative beyond this point. 

In the circular array all density fluctuations move with 
the same sound velocity, so the fluctuation profile does 
not evolve in time and thus overall rotation does not af- 
fect the noise of current / [vg will, however, be important 
later for the analysis of the linear array). To study the 
relaxation of density fluctuations we need to consider the 
deviations of w,. 



5v 



-V5p. 



(20) 



(Including the factor 2 following from Eq. ( |19D would be 
an overestimate of our accuracy, since such nonlinear ve- 
locity can be defined in various ways leading to different 
numerical coefficients). 



where the numerical factor C can be found by comparison 
with the Monte-Carlo results (Figs. ^ and ^), giving a 
value between 1.1 and 1.2. 

Notice that Eq. (|l|) is accurate only if both TV and M 
are sufficiently large. Figure M shows Si{uj) normalized 
by the value 2e(/)(w/27rr)-i'^7V~i[p(l - p)]^/^ for the 
array with TV = 80 and different M. Even at this value of 
TV the plateau corresponding to Eq. ( [2l| ) is not yet very 
wide. With decreasing M the plateau shrinks and there is 
a noticeable deviation from Eq. (HJ). Nevertheless, the 
numerical results presented in Fig. ^ generally confirm 
the analytical result. 

Comparing Eq. ( pl| ) with Eq. ( pj| ) it is simple to esti- 
mate the frequency of the crossover to the high-frequency 
limit: ujh/^TT ^ r[p(l — p)]^, which coincides with the fre- 
quency scale of "collisions" of an electron (or a hole) with 
its neighbors. Notice that for long arrays (TV ^ 1) the 
high-frequency crossover shape does not depend on TV 
(similarly to the linear array case - see Fig. p^ . 

At low frequency Eq. (Efl) becomes invalid when the 
size scale L corresponding to the frequency lul becomes 
comparable with the total array length TV. This al- 
lows us to estimate the position of the low-frequency 
crossover: cj;/27r ~ ^[^(1 - p)]'^/'^N-^/'^, where C is a 
numerical factor. So, we have explained the dependence 
uji cx N~^'^ seen in Fig. |. One can also check that at 
this frequency the result given by Eq. ( |2l] ) transforms 
into Eq. (|l|). 

It is interesting to find out at which electron con- 
centration the single-particle result F = 1/TV becomjes. 
invalid. For TV ^ 1 and small number of electronsEII 
F ~ (Af!)222^^-7(2TV/)!TV, so that considerable devia- 
tion from the single-particle result starts already from 
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M = 2 and scales as M'^/'^. This reflects the fact that in 
ID arrays, significant correlation of hops starts at very 
small concentrations because randomly drifting electrons 
cannot pass each other. 



D. Temperature effect 

In the case of finite temperature when ^ r+, the 
population of all charge configurations remain equal, so 
the average currents satisfy equation 



(/±) = er±p(l~p). 



(22) 



and the net current (/) = (/+) — (/ ) is still given by 
Eq. (H) with r = r+ - r-. Plugging it into Eq. (|), we 
get 



Si{^)/2e{I) = -coth(iy/2T). 



(23) 



This result formally coincides with Eq. (13), but now it 
is only valid for sufficiently high frequencies. 

Figure ||a shows the result of the Monte-Carlo sim- 
ulations for the frequency dependence of the current 
spectral density. As the temperature T is raised be- 
yond the energy difference W , thermal fluctuations grad- 
ually overwhelm the correlation effects, so that the high- 
frequency plateau described by Eq. (|2^) raises and grad- 
ually "floods" regions of lower and lower frequencies. 
[The fact that the low frequency part of the curve is 
less affected by thermal fluctuations can be interpreted 
as follows. Our arguments for Eqs. (|2|) and (|lj) were 
based only on equal distribution of states and Eq. ( |l^ ) 
for the average current, which both remain unchanged 
for arbitrary temperature. So, as long as the tempera- 
ture is small enough so that Eq. (^6|) is still applicable 
for the analysis of fluctuations at the frequency of in- 
terest, the result is virtually unchanged.] Si{lij) may be 
approximately found as the largest of values given by 
Eq. (|2^ ) and the zero temperature result. As soon as 
T > Tc = W[Np{l - p)]^/^, the fluctuations are essen- 
tially thermal at all frequencies, and the Fano factor is 
given by the Nyquist expression 



F 



2T 
IfW' 



(24) 



Notice that as in the single-particle approximation, at 
TV ^ 1 there is a broad temperature region [WN^/"^ ^ 
T <^ WN) where the array is in thermal equilibrium, 
while the Fano factor is still much less than 1. 



E. Coulomb interaction effects 

Coulomb interaction reduces the concentration fluctu- 
ations, so one could also expect a decrease of the current 
fluctuations. This is illustrated in Fig. |6| which shows 



typical Monte-Carlo results for zero temperature. One 
can see that as soon as e^/a becomes comparable or 
larger than W, the low-frequency fluctuations are gradu- 
ally suppressed and can closely approach the limit ([Ti|). 
Figure ^ shows a typical dependence of the Fano factor 
on the array length N for moderate values of the rela- 
tive Coulomb interaction strength r = e^/Wa. At rela- 
tively small TV the scaling F{N) is in between and 
N~^/'^, while eventually at large N it reaches the depen- 
dence F oc N~'^/'^ similar to the case without Coulomb 
interaction. The presence of this transition is specific for 
ID case, since in ID systems the Coulomb interaction 
cannot provide long-range electroneutrality (because the 
electric field {pLe)/L'^ produced by a charged fragment 
of length L, decreases with L). Hence, at large scale 
the density fluctuations are Coulomb-decoupled, which 
makes the general idea of the Fano factor derivation in 
Subsection C valid, leading to the scaling F oc N~^/^. (In 
contrast, in 3D case the Coulomb interaction does pro- 
vide effective long-range electroneutrality, so F inversely 
proportional to the system size is expected.) 

Stronger Coulomb interaction (r > [min(p, 1 — p)]^'^/2) 
tries to fix the distance between the neighboring electrons 
and to turn them into a ID Wigner crystal which may 
be rotated by the external field W. The Fano factor be- 
havior in this case may be rather complex, because it 
depends on whether the integers M and N are "com- 
mensurate" (more strictly, whether their greatest com- 
mon divisor is larger than 1) - see Fig. ^ If it is, beyond 
some critical value Tc of the ratio (about 2.6 for = 20 
and AI = 10, see Fig. ||) the Wigner crystal is stalled (at 
T = 0), the system essentially turning into a Mott dielec- 
tric. At r a little less than Tc the Fano factor starts to 
increase rapidly from F > 1/A^ to some value F^, above 
Tc the ratio F = Si{0)/2e{I) is undetermined, since at 
T = there are neither fiuctuations nor current. In the 
opposite case of "incommensurate" M and N (the g.c.d. 
of M and A is 1) the Mott transition may be absent at 
T — even for arbitrary large r, and both the current 
and Fano factor may tend to the single particle results 
( |l^ ) and ([l3|), respectively. It is curious that on the way 
to this limit the function F(r) may make a bump as if 
it tried to mimic the behavior of its commensurate coun- 
terpart - see Fig. B 



IV. LINEAR ARRAY 

A. The model 

The main change associated with the linear array with 
external electrodes (Fig. |l|a) is that the number M of 
particles in the array is not more fixed. Instead, what is 
fixed are the chemical potentials of the metallic electrodes 
fiL,R relative to the localized state energy. A model of 
the linear array should use this condition to specify rates 
of electron hopping between the electrodes and the edge 
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localized sites. A reasonable way to reduce the number 
of additional parameters is to introduce two extra "edge" 
sites (i = 0, iV, not shown in Fig. |l|a) which are very close 
to the electrodes. Then the "edge" tunneling rates ^ 
are much higher than the "bulk" rates F^, so that the 
edge sites are in thermal equilibrium with the electrodes, 
and the probability of their occupation may be considered 
fractional but fixed: Jl.r = [1 + exp(— yLt^^ij/T)]"-^. In 
this approximation, for a uniform array the rates of tun- 
neling between the edge sites and their neighbors (i = 1 
and i — N — 1) are related to the bulk rates F^ as follows: 



rr = /Lr+, F+ 
rr = (i-/L)r-, 



(i-/i^)r+, 



(25) 
(26) 



We will be interested in the case of identical localized 
sites and similar electrodes, so that — fJ-R and Jl — 
Jr = /■ 

External electrodes also modify the Coulomb interac- 
tion of electrons. Besides that, the image charge effect 
makes the self-energy of the sites dependent on their loca- 
tion, leading to nonuniform transport conditions. Since 
in the present paper we concentrate on uniform arrays, 
we will limit ourselves to the case of negligible Coulomb 
interaction. 



B. Global electron number fluctuation effects 

For the case of T = 0, our modd is reduced to the 
ASEP model with open boundaries. E3 Transport proper- 
ties for the latter model have been studied in detail, es- 
pecially for /l = //J = /• In this case, the probability of 
any charge configuration is the samec3 as if each site had 
independent occupation with probability /. As a conse- 
quence, the dc current is given by Eq. (Hq) with pp./. 
The Fano factor can also be calculated analytically:E3 



= 1-2/(1-/) ^ 



k=0 



k\{k + iy. 



(27) 



and for iV — > oo one finds a simple result:! 

11-2/1, 



(28) 



showing that the shot noise is much higher than in the 
circular arrays - cf. Eq. ([l8|). Only in the evidently spe- 
cial point / = 112, the Fano factor scales as in the closed 
boundary casec3: 



F = (7r7V)-i/2, TV > 1. 



(29) 



Figure ^ (for / = 0.3 and several values of N) shows 
that Si{uj) smoothly decreases with frequency from the 
value given by Eq. (^8|) and eventually reaches the level 
5/(oo) = 2eI/N, in accordance with the general Eq. (|^). 
As we will see later, at large N the frequency dependence 
is quite rich and exhibits three crossovers. 



The fact that the low-frequency shot noise in the lin- 
ear array (Fig. m) at / ^ 1/2 is much higher than in the 
ring array (FigTllp) has a simple explanation: the total 
number of electrons in the case of "open boundary condi- 
tions" may significantly fluctuate, while on the ring this 
number is fixed. Analytically, this effect may be espe- 
cially simply considered for the case / ^ 1 (or similarly 
1-/^1). Then the array is empty most of the time, 
and is entered very rarely by an electron (or hole). Af- 
ter the entry, the electron is transferred in a succession 
of hops through the array, in total transferring charge e 
from one electrode to another. This is exactly the situa- 
tion for which the original Schottky formula was derived, 
so that we get F = 1 in agreement with the correspond- 
ing limit of Eq. (H). 

In the case / <C N^^^"^ (when electrons do not col- 
lide with each other) the frequencv-dependence of Si can 
be obtained from Eq. (14) of RefEj, which was derived 
from the orthodox theory of single-electron tunneling for 
the similar sequential transport scenario. Assuming that 
tunneling rates are equal, F^ = F, besides the negligibly 
small rate F]*", we get 



2el 



1 

N 



7p 



Re(l -iw/F) 



(30) 



we have confirmed this result using Monte-Carlo simula- 
tions. For ^ 1 this formula is reduced to 



Siju;) 
2el 



sin(Afw/2F) 
Nu;/2T 



(31) 



that is obviously the normalized and squared Fourier im- 
age of the rectangular envelope of the train of N current 
pulses during the single electron passage. 

For / ^ 1, Eq. ( p8| ) may be interpreted as follows. 
Let us again apply Eq. (|l6|) to long-range fluctuations. 
Then since dl/dp = eF(l - 2p), one finds 5/(0) = 
g2p2j|2 _ 2p)^Sp{0), where S'p(O) is the low-frequency 
intensity of fluctuations of the total array occupation 
(p EE M/N). Notice that for p = / = 1/2 the result 
vanishes, and we should go after the higher order effect 
as we did for the ring array. For all other values of /, 
we may use the estimate S'p(O) {{S p)"^) / Auj , where 
(((5p)2) = /(I _ /)/iV and the effective bandwidth Acu 
can be estimated as |ds/A^| (unlike in the ring array, the 
density fluctuation is carried out of the linear array with 
velocity Vs given by Eq. (|l^)). Combining these formulas, 
we obtain Aw - 11-2/1 F/7V,5p(0)^/(l-/)/F|l-2/|, 
and F — const x |1 — 2/|. The numerical factor in this 
result for the Fano factor cannot be derived in this crude 
way, but it obviously equals unity because at / = 1 we 
should get the previous result, F — 1. Thus we com- 
pletely recover the exact result (p8|). 

One more possible derivation of that equation can be 
obtained along the following line. If / < 1/2, then the 
electrons can be supplied from the left electrode with 
the maximum rate /F, while the average "sink" velocity 
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(1 — /)r is larger. Hence, only electrons relatively close 
to the left boundary can affect the entrance of the next 
electrons, and so the low-frequency correlation is essen- 
tially the boundary effect. Using this idea and taking into 
account, for example, correlations only due to the three 
first jumps, it is easy to obtain F = 1 — 2f + 0{f^). Tak- 
ing into account more jumps we would eventually show 
that Eq. ( p8| ) is exact. 

A concentration fluctuation supplied from the bound- 
ary moves with velocity Vs, so for ^ 1 the correspond- 
ing envelope in I{t) has rectangular shape with duration 
7V/|ws|. Combining the corresponding frequency depen- 
dence of Si{uj) with the exact result for F, we get 



2el 



11-2/1 



sin(A^cj/2r|l - 2/1) 



Nuj/2T\l - 2/1 



(32) 



At sufficiently high frequency, cu 3> \vs\/N, the "non- 
linear" contribution from the concentration fluctuations 
obviously should be the same in the linear and ring arrays 
(with equal average concentration p = f). Hence, Si{ll!) 
will still be given by Eq. (|2^) while in the crossover region 
it can be crudely estimated as a sum (or maximum value) 
of two contributions given by Eqs. (g|) and (^. As a 
result, there are three characteristic frequencies in Si{lj) 
dependence at N ^ 1: the low- frequency saturation oc- 
curs at w < ~ r|l — 2/l/Af, the intermediate- frequency 
dependence described by Eq. ( ^ starts at uj > uj„i ^ 
rAr-3/5|i _ 2/|9/5[/(l - /)]-2/5, and finally the high fre- 
quency saturation occurs at uj ^ Uh ^ r[/(l — /)]'^> sim- 
ilar to the ring array case. 

The case / = 1/2 plays a special role in the ASEP 
theory, as can be easily noticed comparing Eqs. ( p^ ) and 
(p9|). Actually, this case is quite important since for suf- 
ficiently long arrays with Jl > 1/2 and Jr < 1/2 tjia 
electron concentration in the bulk of the array is closeEII 
to / — 1/2 (so (/} = er/4) and, hence, the scaling 
F oc N^^/"^ holds as in Eg-, (p9|) . (As an example, for 
/l = 1, /i? = the result IsEI = 3(27r)i/Vl67Vi/2). At 
f — 1/2 the low-frequency fluctuations can no longer be 
considered as a boundary effect, since the "sink" velocity 
(1 — /)r is equal in this case to the maximum supply 
rate /F; hence, the transport becomes jammed and the 
correlations involve the whole array length. 

In the respect that the boundary effects are no longer 
important, the linear array at / = 1/2 is very similar to 
the circle array. Figure ^ shows the frequency depen- 
dence of the current spectral density for / = 1/2 and 
several values of N . The data look similar to that in Fig. 
^. The main feature is dependence in the interme- 

diate frequency range. To check the validity of E q. (Ell) in 
this range. Fig. |l0| b shows the same data as Fig. |l0|a"but 
normalized by Sq{uj) = 2el{uj /2'K)~^I^N^^[f{l - f)fl^. 
We see that as N grows, the intermediate region becomes 
more and more pronounced. 



C. Temperature effects 

Figure ^ shows the numerically calculated effect of 
nonvanishing temperature on the shot noise in a linear 
array. It shows that the effect is quite similar to that 
in a ring array (Fig. ||), however, because of the higher 
initial intensity of low-frequency fluctuations (at T — Q) 
the noise becomes completely thermal at a higher tem- 
perature, T > VF|1 - 2/|. 



V. DISCUSSION 

Probably the most important result of our analysis is 
that in contrast to the expectation based on the analysis 
of ID arrays of conventional tunnel junctions, the shot 
noise in uniform ID hopping arrays is typically much 
higher than 1/A^ of the Schottky value Si — 2e{I). How- 
ever, in some cases this lower bound can be achieved. In 
order to sort out these cases, it is useful to consider the 
current fluctuations as the result of three major sources: 

- time randomness of electron tunneling events, 

- electron density fluctuations, and 

- thermal fluctuations. 

Crudely speaking, the lower bound 2e{I)/N for the noise 
is determined by the first contribution, while the second 
contribution typically increases the noise significantly 
even at T = 0. 

At relatively high frequencies the current spectral den- 
sity in a ring array and ID array between electrodes be- 
haves pretty similarly. In particular, the high-frequency 
asymptote is given by the same Eq. (|^) and is determined 
by capacitive factors A,;. (This result is also valid for the 
conventional case of ID array of tunnel junctions - see 
Appendix). If Xi = 1/iV, then in all cases at low temper- 
ature T we have Si{oo) = 2e{I)/N. 

However, at low frequency the noise behavior in a ring 
array and a linear array is quite different. The reason 
for the difference is that in a linear array the total num- 
ber of electrons can fluctuate while in the ring array it 
is fixed. In the case when the single-particle approxima- 
tion is applicable for a linear array, the relative density 
fluctuations are maximal, and the Fano factor is not sup- 
pressed: F = 1 at T = 0. The electron "collisions" (the 
Pauli exclusion) reduce these fluctuations, but quite inef- 
ficiently. Only in the special case of half filling (/ = 0.5) 
when "traffic jams" have all size scales, the Fano factor 
decreases as N"^/-^ with the array length N; in other 
cases the dependence F{N) quickly saturates at the level 
F = |1 — 2/|. One can speculate that Coulomb interac- 
tion should be a more efficient factor in suppression of F, 
since it may significantly reduce the electron density fiuc- 
tuations, however, this conclusion has still to be verified 
numerically. 

In contrast to the linear array, in the uniform ring ar- 
ray the uncorrelated motion (of a single electron) pro- 
vides the maximal suppression of the Fano factor, F = 
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The Pauli exclusion in fact increases F leading 
to F cx N~^^^. However, the extra correlations due to 
Coulomb interaction between electrons on different sites 
make transport "smoother" and reduce the Fano factor, 
in some cases down to the lower bound F — 1/N. 

It is instructive to compare these results with those for 
a ID array of tunnel junctions (see Appendix). In the 
latter model the Fano factor is determined purely by the 
junction resistances. In some sense, this is a consequence 
of strong Coulomb interaction which forbids noticeable 
charge fluctuations and establishes fast long-range corre- 
lations between currents through different tunnel junc- 
tions. In the uniform array at low temperature the noise 
suppression is maximal, F — 1/N. However, if the junc- 
tions are very small, single-electron effects can lead to 
significant charge fluctuations and thus, increase the Fano 
factor. For example, F = 1 is realizedo in the vicinity of 
the Coulomb blockade threshold when the transport has 
a bottleneck even in the uniform array. 

So far we have reviewed our results for the uniform 
case. Now let us briefly discuss hopping transport noise 
in nonuniform ID arrays. It is simple to study one par- 
ticle inside the ring array with arbitrary disorder at low 
temperatures. In this case the transport is unidirectional, 
= 0, and the average current is obviously given by the 
expression 



and 1. The minimal rate ("bottleneck") Tmin will be 
about ro(2ao/^eAf)2°o/« on average (here e — 2.71..), 
while the next minimal rate Tmin+i will be much larger, 
r,nin+i/Trmn ~ (2ao/^eiV)-^''o/« > 1. It is casy to see 
that in this case both the average current [Eq. (^)] and 
the Fano factor [Eq. (^5|)] are determined by the bottle- 
neck: (/) = eVmin and F = 1. 

It is also instructive to consider a model where the 
maximal distance is limited by some big value amax 
{amax ^ clq)- For example, this describes the situation 
in which some other transport mechanism starts to dom- 
inate over tunneling when the sites are too far apart, 
thus limiting F from below. If TV ^ exp(amaa;/ao)j 
the results for the average current and the Fano factor 
do not differ from the case considered above. However, 
for very long arrays, N ^ exp(omaa:/ao)7 the transport 
is limited by many similar bottlenecks, so that (/) ~ 
eFo eyi'p{—2amaxl^-'rO-max/oLVi)/N and the Fano factor de- 
creases with the array length, F ~ e'x.p{amax / ao) / N . 

We can use this result for a preliminary estimate of the 
Fano factor for hopping in disordered 2D and 3D systems 
when the-jtransport is mainly determined by percolation 
clusters.Ell If the single-particle approach is applicable, 
then as in the case above we can simply count the number 
of similar bottlenecks in the transport direction. With 
this argumentation, we obtain a simple estimate 



{I)=e 



(33) 



while the formula ijes the current spectral density has 
been derived in Rcf.E3 for the case Xi — 1/N, and can be 
readily generalized to include arbitrary A,;: 



N 



5/(w) = 2e(/) ^ -I- 4:6(1} Re 



N 



1=1 

N-1 N 



N 



Ll=l 



1 + 



.1=1 



1 = 1 m=l 



k=l 



r+ 

^ k+m 



, (34) 



where by definition F^^ 
formula is reduced to 



F^. At zero frequency this 



5/(0) 

MI) 



(35) 



F^LJL, 



(36) 



where L is the sample length and Lc is the characteristic 
size of the percolation cluster. The applicability range of 
the approach is rather unclear, so this result still has to 
be confirmed using either more quantitative analysis or 
numerical Monte-Carlo modeling. (Preliminary numer- 
ical analysis of hopping in uniform 2D arrays gives an 
indication of the dependence F oc L~'^ where v ^ 0.8 is 
rather close to unity.) 

In conclusion, in the present paper we have studied 
the shot noise at hopping in ID arrays of sites, concen- 
trating on the uniform case and briefiy considering the 
effect of disorder. It is important to extend this study 
to 2D and 3D hopping. The presented results hint that 
the Coulomb interaction may play the crucial role in the 
suppression of low-frequency shot noise. 
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and allows study of the statistics of the Fano factor for 
random distribution of F^ in a long array, ^ 1. 

As the major factor, let us take into account the de- 
pendence of tunneling rate on the distance between 
sites, F^ = Fq exp(— 2ai/^), where ^ is the localization 
length, and assume that independent random obey 



the Poisson distribution, p{ai 



exp(— ai/ao) where 



oq ^ ^ is the average spacing. Then the distribution 
of rates can be parameterized as F^ — T^xf"^^^, where 
the random number Xj has uniform distribution between 
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APPENDIX A: ID ARRAY OF TUNNEL 
JUNCTIONS 

In this Appendix we calculate the spectral density of 
the shot noise in a ID array of tunnel junctions. We 
will first use the standard circuit theory and then ex- 
tend the calculations to the case of weak single-electron 
effects (which are assumed to be small because of high 
temperature or high current). 

Figure [ij shows the array of N tunnel junctions in se- 
ries. Each junction is characterized by resistance Ri and 
capacitance d {i = 1, . . . N), while Cg^k {k = 1, . . . N — 1) 
denote capacitances of "islands" to the ground. This cir- 
cuit does not describe the long-range capacitances which 
can be espiecially important for islands and junctions of 
small size,E3 that can require, aumerical calculation of 
the total capacitance matrixEjES The derivation below 
is valid for arbitrary capacitance matrix, however, for 
simplicity we will refer to Fig. ^ 

In the absence of single-electron correlations the I-V 
curve of the array is hnear, (/) = V/R^, i?s — J2i-^ij 
and the average voltage across each junction is propor- 
tional to the junction resistance, (Vi) — {I)Ri. To stucb, 
the fluctuations we follow the standard circuit theorjO 
and introduce the sources of the current noise (t) in par- 
allel with the junctions. Separating the current through 
ith junction into two parts flowing in opposite directions, 

{i) = {i^)-{in, (i+)/{in-^M-e{v.)m, (ai) 

and using the Schottky formula for each part, we get the 
following white spectral density for the noise source (t) : 

S^i{uj) = 2e(/+) + 2e{ir) = 2e(/) coth{e{Vi) /2T). (A2) 

Let us denote by (piit) the fluctuating part of the ith 
island potential {(j)o{t) — (pNit) — because we assume 
constant potentials of the leads), then the current Ii{t) 
through ith junction can be written as 

= (/) + [(b,-i{t) - Mt)]/R^ + ^^it)■ (A3) 
The evolution of 4>i{t) is described by the equation 

N 



(A4) 



where D = C^^ is the inverse capacitance matrix [ob- 
viously Di^j = Dj,i and i'o.i = Dna — 0] and can be 
rewritten in the following form: 



JV-l 



N 



k=l k=l 
Ai,k = Bi^k+l/Rk+1 — Bi^k/Rk, 

Bi^k = A,/c — A,fc-i- 



(A5) 

(A6) 
(A7) 



Notice that A is (iV - 1) x (TV - 1) matrix while B is 
{N — 1) X N matrix. In the frequency representation 



Eq. ( [A5| ) can be written as lui 4>{^) = A 0(w) -I- B ^(w) (i 
is the imaginary unity) and can be easily solved in the 
matrix form. 



</>(c.) 



;1 - A 



B^(c.). 



(A8) 



Using Eq. (A2) we find the Fourier transform of liit): 



N 



(A9) 



N-l 



k=l L 



- A 



i-l.k 



A 



B 



(AlO) 



where by definition (ilo — A] 



0,k 



(tuj - A 



0. 



No tice that at w — > oo the only surviving term in Eq. 
( AlC ) is the Kronecker symbol Sij so that 



(All) 



At a; = it is possible to prove (using somewhat cum- 
bersome algebra) the relation 



/.(0) = ^(i?,/i?E)^,(0), 



(A12) 



which obviously means that at low frequencies the cur- 
rent is distributed accor ding to resistances and equal in 
all junctions. (Eq. ( A12 ) shows that the fraction Rj/R^ 
of the current flows through the array while the rest 
is returned via the "shunt" Rj). 

The spectral density of the current /; through ith junc- 
tion can be readily calculated as 



(A13) 



because the nois e s ourc es S,i ar e mu tually uncorrelated. 
Using Eqs. (^), ( |All| ) and, ( |A12| ) we get the simple 
expressions in the limiting cases: 

Su{0) = 2e{I) ^(i?2/i?|) coth(e(/)i?,/2r), (A14) 



SI^{oo) = 2e{I) coth(e(/)i?i/2T). 



(A15) 



In experiment it is usually impossible to measure the 
current through one junction, and the only measurable 
quantity is the current in the external lead which con- 
tains the contribution from the displacement current and, 
hence, depends on the currents through all junctions. 
The current I{t) at the left external lead can be expressed 
as the linear combination, 

/(<) = ^A,/,(<), (A16) 

i 

K^5u + Y,Co,jB,^,, Y,X, = l, (A17) 
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where Co,fe is the element (always negative) of the ca- 
pacitance matrix between the left electrode and fcth is- 
land. In the case of Fig. |l2|, Co,fe = —CiSik because 
the left electrode is capacitively coupled only with the 
first island, so \i — — CiBi i 



If aU Cg,, = in 
Fig. then — C^^/J2i ^7^^ while for a long uni- 



form array, Ci = C = const, Cg,i = Cg = const, 
> C/Cg, the coefficients are = X'-'^{l-X) where 
X = 1 + Cg/2C - [(Cg/2C)2 + If thc tunnel- 

ing system is different from what is shown in Fig. ^ and 
consists of small islands which are separated by much 
larger distances and imbedded into the plane capacitor 
(external electrodes), then similar to the hopping case 
Ai — ciil "Ylij o-j where Ui is the length of the projection of 
ith tunneling jump onto the direction perpendicular to 
the capacitor planes. 

Using Eqs. (Al(;) and ( [A9| ) we obtain the following 
expression for the spectral density of the left external 
current: 



5/M = E 



s> 



(A18) 



(notice that at finite frequency, 5/(a;) for the left and 
right electrodes can be different). In the particular case 
shown in Fig. n2. 



N-l r 



i 771 , k 



— h Ci^i^rn 

Hi 



(A19) 



In the limit a; = the spectral density >S'7(0) of the 
external current coincides with the spectral density of 
the current through an y jun ction (see Section II) and, 
hence, is given by Eq. (A14). While at low frequency 
Si is determined only by the circuit resistances, in the 
opposite limit, cj ^ oo, it is determined mainly by the 
circuit capacitances, 

5/(oo) =Y,^^t % = 2e(/) coth(e(/)i?,;/2T) 

i i 

(A20) 

(resistances are important only when the Nyquist noise 
contribution is considerable, IRi ^ T). It is sim- 
ple to check that in the thermal equ ilibrium, / = 0, 
the low- frequency noise [sec Eq. (A14)] always satisfies 
the Nyquist formula, 5/(0) — AT / R^, however, this 
is not true at finite frequency, for example 5/(oo) = 
4T^^ X^/Ri. Notice that in this formalism u = oo still 
means huj <C max(T, e^i). However, it would be simple 
to take into account zero-point fiuctuations by replacing 
Eq. (^) withka %(u;) = R'^ E±(e^^ ± ^w) coth[(eV - ± 
noj)/2T\. 

In a uniform array at zero temperature the nois e at 
low frequency is suppressed N times [see Eq. (A14)] in 



comparison with the Schottky formula, 5/(0) = 5/^(0) = 
2eI/N. However, at high frequency the suppression of 
the external current noise is usually weaker, 5/(oo) = 
2e/E»\^ < EjA? < 1> while for the current 

through a particular junction there is no suppression at 
all, Sn{oo) = 2el. 

It is interesting to find out at which lo the low- 
frequency result is no longer accurate. We have stud- 
ied the long uniform arrays shown in Fig. N ^ 
[C/CgY/"^, and found numerically that the relative ac- 
curacy e <C 1 of the low-frequency result, 5/(u;)/2e/ = 
(1 + corresponds to the frequency uj/2it ~ 

l.l{RCg)^^e^/^N^^. This formula, however, cannot be 
used to describe the crossover to the high-frequency 
asymptote, 5/(oo) = (1 + AC /Cg 



^1/2 



If the typical sizCpOf the tunnel junctions is small, 
single-electron effectfj become important. Below we 
extend the standard noise theory to this case, assum- 
ing that single-electron effects are weak due to relatively 
high temperature, T > e^/C, or relatively high current, 
/ > e/RC {C and R ^ Rq are the typical capacitance 
and resistance of tunnel junctions). j— . 

According to the "orthodox" theory,Ej the rate of tun- 
neling through ith junction (in the positive direction). 



r = V^^^/eR^l - exp(-ey/-^VT) 



(A21) 



is governed by the effective voltage V^'^'^ which is always 
smaller than the actual voltage Vi, 



= V,- e/2Ct,r, 



(A22) 
(A23) 



where Ct^i is the total capacitance of ith junction. (For 
tunneling in the opposite direction the effective voltage 
is —Vi — e/2Ct,i-) Linearizing Eq. (A21) and averaging 



over the fluctuating Vi{t), we obtain the equation 

(!) = - {in, 

±{V) - e/2Ct,^ 



(It) = 



R,[l - exp(-e(±(l/,) - e/2Ct,^)/T)] ' 



(A24) 
(A25) 



which allows the calculation of average voltages (V^) 
for the given average current (/). (Notice that the 
high- voltage offset of the /- V curve is equal to Vof / = 
J2i^/'^^t,i exactly.) This approximation has been suc- 
cessfully used for the analytical and numerical analy|Sis 
of the /-V curves of Coulomb blockade thermometers.E3 
Since we assumed an essentially linear response of any 
junction current, it is natural to use the formalism of 
the standard circuit theory, so the result for the cur- 
rent spectral densit y wi ll be g iven by the same Eqs. 
(|aTc| ), ( [aT^ ), and (|Al7D-(|Al9f). However, the second 
equality in Eq. ( |A2| ) is no longer valid, and we have at 
least three choices for the "seed" noise spectral density: 
= 2e(/+) + 2e(/-), S^^ = 2e(/)coth(e(y,)/2r), or 
5{i = 2e(J)coth(e(J)i?,/2r). The flrst choice seems to 
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be the most natural one, however, numerical comparison 
with the results of Monte-Carlo simulations shows that 
the first formula usually underestimates noise (see Fig. 
[l3| ), the third formula overestimates it, and the second 
formula (which is in between two others) usually gives the 
closer result, though not always. Notice, however, that 
all three approximations coincide in the limits of both low 
and high temperature, so the difference between them is 
never too significant within the applicability range of the 
formalism. 

As the next level of approximation for the "seed" noise 
Sf^i at finite temperature, it is possible to estimate the 
standard deviation of fluctuating Vi (t) and take into ac- 
count the correction due to the second derivative of Eq. 
( A-21 ) . It is also possible to take into account the effective 
increase o f th e junction resistances used in the evolution 
equation ( |A3| ) at finite temperatures (we did not imple- 
ment this last idea numerically). Figure ^ shows the fre- 
quency dependence of the spectral density of the current 
in the (left) external electrode for the uniform array of 
iV = 10 junctions with Cg/C = 0.1 symmetrically biased 
by voltage V = 5e/C at temperature T = e^/C. The 
thick line shows the results of Monte-Carlo simulations 
while t he th in lines represent the calculations using Eqs. 
( Al^ )-( A19). For the lowest and highest thin lines, the 
first and third formulas for S^i discussed above have been 
used. The thin line corresponding to the second formula 
is almost indistinguishable from the thin line showing the 
result using the second (quadratic) approximation of 5*^^. 
As one can see, these lines are quite close to the Monte- 
Carlo result. With an increase of current or temperature 
all thin lines become closer to each other, and the agree- 
ment with the Monte-Carlo result becomes even better. 
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FIG. 1. (a) Linear array of N ~1 localized sites connecting 
two electrodes ("open boundary conditions"). The electron 
transport is determined by the tunneling rates F^. (b) Cir- 
cular array ("periodic boundary conditions") with A*' sites 
occupied by M electrons. 

FIG. 2. The Fano factor _F as a function of energy differ- 
ence per site in a circular array occupied by one electron. 

FIG. 3. Frequency dependence of the spectral density 
Si{ij) for uniform circular arrays at T = for several val- 
ues of array lengths TV and electron concentration p = M/N: 
(a) p = 0.3, (b) p = 0.5. 

FIG. 4. Current spectral density for the array in a circle 
normalized by Sq{uj) = 2eI{Lo N-^[p{l - p)f^^ [see 
Eq. 

FIG. 5. Current spectral density for the ring array with 
A'' = 80 and M = 40 for several temperatures. 

FIG. 6. Current spectral density in a ring array at 
zero temperature for several values of Coulomb interaction 
strength r = /aW . 

FIG. 7. Fano factor of ring arrays with a fixed electron 
concentration (M/N — 0.5) as function of the array length 
A^, for several values of Coulomb interaction strength. Lines 
are just guides for the eye. 

FIG. 8. The dependence of the Fano factor for the ring 
array on the strength r of Coulomb interaction. 



FIG. 9. Frequency dependence of the spectral density 
Si{ij) for uniform linear arrays with symmetric boundary con- 
ditions, fL = fR = 0.3, at T = 0. 



FIG. 10. (a) - Current spectral density for linear array 
with /l ~ fa ~ 0.5 at T = 0. Notice the dependence 
Si{u)) ~ Lu"^^'"^ in the intermediate frequency range between 
the saturations at low frequency {F ~ N~^^^) and high fre- 
quency {Si{oo) /2el = 1/A'^). (b) - The same data normalized 
by SoH = 2eI{c,/27v)-'^'N-'[f{l ~ j)fl\ 

FIG. 11. Frequency dependence of the current spectral 
density in the linear array with A*' = 10 and /i, = fn — 0.3 at 
several temparatures. 

FIG. 12. ID array of N tunnel junctions with capacitances 
Ci and resistances Ri. Cg^i are the capacitances to the ground. 
The Langevin noise sources are presented by random current 
generators £,i{t) parallel to the junctions. 



FIG. 13. Normalized spectral density Si{u;) for a uniform 
array of small tunnel junctions. Thick line shows the result of 
Monte-Carl o sim ul ation s while thin lines 1-4 are calculated 
using Eqs. ( |A18[ )-( |aT£| ). For Unes 1-3 the "seed" noise 5"^^ 
is calculated as 2e{{lf) - (Ir)), 2e(/) coth(e(K)/2T), and 
2e{I) coth{e{I)Ri/2T), respectively. For the line 4 (by coinci- 
dence, almost indistinguishable from line 3) S^i is calculated 
using the quadratic approximation. 
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